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Abstract 

We treat the problem of N interacting, axisymmetric black holes 
and obtain two relations among physical parameters of the system in- 
cluding the force between the black holes. The first relation involves 
the total mass, the angular momenta, the distances and the forces be- 
tween the black holes. The second one relates the angular momentum 
and area of each black hole with the forces acting on it. 

1 Introduction 

The problem of interacting black holes dates back to the beginnings of Gen- 
eral Relativity. Although a solution representing multiple black holes in 
equilibrium is known to exist [3] (the Majumdar-Papapetrou solution, which 
consists of N extreme Reissner-Nordstrom black holes, i.e. each black hole 
has electric charge equal to its mass parameter, = m,), it is in general ex- 
pected that there should be a force between different components of the hori- 
zon in order to prevent the spacetime from collapsing. Nevertheless, based 
on Newtonian ideas, a main concern was the following: Given two black 
holes, could the spin-spin repulsion compensate for the gravitational attrac- 
tion and keep the system in equilibrium? Many attempts have since then 
been made in order to answer this question (see [2], |28| . |18| |24] . j9], |21]). 
In particular, Li and Tian [18] and Neugebauer and Hennig [22j, [16j |23| 
proved that an axially symmetric, stationary, vacuum solution of Einstein 
equation with disconnected horizon (and only two connected components) 
can not exist in equilibrium. Weinstein [28], [26], |27J, using harmonic maps, 
showed that asymptotically flat, vacuum, axially symmetric multi black hole 
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solutions exist, although possibly with a conical singularity in the bounded 
component of the symmetry axis. 

The appearence of this conical singularity on the axis was clearly seen 
in the superposition of two Schwarzschild black holes in the early work of 
Bach [2] and was later discussed by Einstein and Rosen in [llj. Although 
the existence of a singularity may seem discouraging, the solution so built 
has a well defined action and the corresponding thermodynamical potential 
describes a system of black holes at fixed temperature and at a fixed distance 
|15| . |6]. The long range interactions in the multiple black hole solutions are 
unable to provide equilibrium between the various black holes and the conical 
singularity is interpreted as a boundary condition that keep the holes at fixed 
separation preventing the collapse of the system. 

The non-existence proof of Neugebauer and Hennig [22] for two rotating 
black holes in vacuum strongly uses the relation A > 8ir\J\ between the area 
and the angular momentum of each black hole which is known to hold in the 
axisymmetric, regular case (see pQ). This type of geometrical inequalities 
relating physical parameters of black holes have been extensively studied in 
the last few years as an attempt to explore the allowed values of angular 
momentum, mass and size a black hole can have, and turned out to be 
relevant in different contexts as well (see [7] for an excellent recent review on 
geometrical inequalities). As interesting outcomes of these researches, two 
main relations have been found. 

The first one of such relations in the context of multiple black holes con- 
figurations was proved by Chrusciel et al [5] (see also |10] for the inequality 
in the case of AT = 1), and gives a lower bound to the ADM mass in the 
axisymmetric, vacuum case in the form 



where m is the total (ADM) mass, Jj are the Komar angular momenta of the 
individual constituents and / is a function of the J, (and possibly depending 
on other parameters as well), not known explicitely, which reduces in the 
case of a single black hole to v/| J|. 




The second inequality we are interested in, and the one used in [22] is the 
relation 



between the angular momentum J and the area A of an axisymmetric black 
hole which has recently been proven. In [T], [12] and [8] maximal vacuum 
initial data (possibly with a positive cosmological constant) containing a min- 
imal stable surface are considered. Also, in [T7] the inequality is proven for 
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spacetimes having outermost stably marginally trapped surfaces, not neces- 
sarily in vacuum. As we will discuss later, ([2]) has been extended to incorpo- 
rate electric and magnetic charges. Inequality (|2J) holds for each black hole in 
a regular initial data and, loosely speaking, states that a hole can not rotate 
too fast. 

An important ingredient in the above results, relations ([I]) and (J2J), is the 
assumed regularity at the symmetry axis. As we mentioned, this is related 
to the notion of equilibrium. We will incorporate the interaction between the 
black holes into the problem is by means of relaxing the regularity condition 
at the axis and therefore allowing a conical singularity between the holes. 
The value of the deficiency angle will be related to the gravitational force 
|28| . The purpose of this work is to investigate how inequalities ([I]) and (j2J) 
are modified when one takes into account explicitely the interaction between 
the black holes, that is, when one includes non-equilibrium configurations. 
In this way we will obtain bounds on the forces between the black holes (i.e. 
on the deficit angles) set by the physical quantities A, m, and J. 

The article is organized as follows. In section [5] we briefly introduce the 
basic features of multiple black hole solutions to Einstein constraint equations 
and the explicit relation between force and deficit angle. Then, in section 
12.11 we present the more especific hypothesis needed in order to obtain a 
relation between mass, angular momentum and force, similar in espirit to 
([T]) , theorem 12.11 and finally, in section 12.21 we find a new inequality relating 
area, angular momentum and force, theorem 12. 2 \ analogous to (|2J). 

2 Main result 

We begin this section with a brief description of multiple black hole solutions 
of Einstein equations in the form of maximal axisymmetric initial data. 

Consider a maximal, simply connected initial data set (M, g a bK a b), in- 
variant under a U(l) action. It is known [4j that there exists a global coor- 
dinate system (p, z, <fi) such that the metric takes the form 

g = e a+2q (dp 2 + dz 2 ) + pV(# + P A p dp + A z dz) 2 , (3) 

where the functions are taken to be 0— independent. 

We assume vacuum for simplicity, although matter satisfying some posi- 
tive energy condition would be also allowed. The constraint equations in the 
maximal case, g ab K a b = 0, are 

3 R = K ah K a \ W a K ah = 0. (4) 
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Due to these constraints, there exists a twist potential u related to the second 
fundamental form by the expression (see [10]) 



e abc K h e rfrfdx a = du (5) 

where r] a is a 2ir— periodic Killing vector field and e abc is the Levi Civita 
tensor. Moreover, we have the useful bound [10] 

K ab K« b > e— 2 ^. (6) 

The quantities that we are interested in in this work are the total mass 
of the initial data, and the area and angular momentum of each black hole. 
As is well known, the mass is a global quantity, given by a boundary integral 
at spatial infinity 

m = -z- lim / (d b g ab - d a g bb )n a dS (7) 
16vr r^oo J dBr 



where r = i/p 2 + z 2 is the Euclidean distance, d denotes partial derivatives, 
B r is the Euclidean ball r = const, dB r is its border, n a is its exterior unit 
normal and dS is the unit surface element with respect to the Euclidean 
metric, dS = sin OdOdcf). 

Now, consider an arbitrary 2-surface E in M with induced metric 7 a b, 
then the area of E is given by the expression 

A(E) = / dS^ (8) 
is 

with dSry the surface element with respect to 7^. 

Finally, for axially symmetric and maximal data one can define the angu- 
lar momentum J associated with E (the Komar integral of the Killing vector) 
as 

J(E) = J^K abV a n b dS 7 , (9) 

where n a is the unit normal vector to E. 

For concreteness, in what follows, we consider N > 1 black holes located 



at the symmetry axis, which will be represented by punctures in section 12.11 
and by minimal stable surfaces in section 12.21 The axis p = is denoted by 
T. It contains N — 1 bounded components denoted by Tj, % = 1,..,N — 1, 
and two unbounded components, T and r^v- 
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As we mentioned in the introduction, the interaction between the black 
holes can be realised through the non regular character of the solution at T. 
The regularity condition at T can be stated as 



lim 



1 



(10) 



e a / 2 +idp 



where r] = r] a r] aj and is translated into the requirement that the metric func- 
tion q (see the expression ([3])) vanishes at the axis: q\r = 0. 

In general this condition will not be satisfied and the metric will not 
define a regular solution of Einstein equations. If (jlOp does not hold, then 
the singularity at the axis is called a conical singularity The deficit angle at 
the axis is given by e~ Qi where qi is the constant value of the function q at the 
i— th connected component of the axis. These singularities are interpreted as 
forces needed to balance the gravitational attraction and keep the bodies in 
equilibrium and can be expressed in terms of as (see |27] for more details) 



One would expect this force to be positive, reflecting the fact that a positive 
force is needed in order to prevent the black holes from falling on each other. 
A positive force has been found in a few situations, like in the zero |2] and 
small angular momentum case [19], the case when the problem admits an 
involutive symmetry |18j . when there are two equal, counter-rotating Kerr 
black holes [25J and in the limit when one of the black holes becomes ex- 
treme and the distance to one of the adjacent black holes tends to zero [27J. 
Nevertheless, so far there is no general proof stating that the force is always 
positive. 

One can prescribe the value of the function q on one of the component 
of the axis [27J, the other values will be determined by the solution of the 
constraints, and are not known in general. In this article, we give two results, 
which are the extensions of inequalities (pQ) and ([2]) to include the non zero 
values of q at the axis and are presented in the following sections. 

2.1 Bound in terms of mass and angular momentum 

In order to extend Chrusciel et al result [5J, inequality ([T]), to the non-regular 
case, we state here their main assumptions. We will retain basically all of 
them, and only relax the regularity of the function q at the axis. 

As we mentioned before, for the mass-angular momentum inequality we 
represent the black holes as iV punctures on the axis, separated by distances 
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Lj, each puncture being an asymptotically flat end in the sense that for k > 6, 

gab-5 ab = o k (r- 1 / 2 ), (12) 

d k9ab e L 2 (M) (13) 
where / = o k (r x ) means that / satisfies 

d kl ...d kn f = o(r x ~ n ), 0<n<k. (14) 

The main result of this section gives a relation among the total mass of the 
initial data, the individual angular momenta of the each asymptotic region 
and the forces between them. 

Theorem 2.1. Consider a maximal, asymptotically flat, axisymmetric initial 
data (M, g a b, K a b) as described above, with N extra asymptotically flat ends 
represented by punctures on the axis T, with angular momenta J i; i = 1, .., N. 
Then, there is a lower bound to the ADM mass m given by 

N-l 

-J h(l + >f{J x ,...,J N )-m (15) 

-i=i 

where Li is the Newtonian length of Ti, T% is the force between adjacent 
punctures and f is a function of the 's, . 

Before entering into the proof, we want to make a few remarks about 
the result. It is clearly seen from inequality (I15p that it does not imply, in 
general, the result > / — m. That is only the case when the forces are non 
positive. But if the forces are positive, which is the expected situation, in 
principle, the total mass m could take values lower than /, as opposed to the 
regular case [5]. 

As the data are asymptotically flat, we must have go = Qn = 0, and (I15j) 
reduces to Chrusciel et al result, m > f, for = and also to Dain's result, 
m > y/\T\, for iV = 1 [10J. 

Note that there is no rigidity statement in this theorem. This is due in 
part, to the lack of a known minimizer for a mass functional appearing in 
the proof of Chrusciel et al result in the case N > 2. But also, because 
in inequality (|T5|) there appear two quantities that depend on the given pa- 
rameters of the data. Both, the mass and the force are to be computed 
after solving the constraints and can not be prescribed a priori. Therefore, 
for given values of angular momenta Jj and separation distances Lj, there 
might exist different solutions saturating the inequality with different values 
of Ti and mass. In this respect, it is worth mentioning, that there exists an 
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explicit solution describing two equal, counter-rotating extreme Kerr black 
holes given by Manko et al |2T] which could saturate ( IT5|) . We will meet this 
solution again in the following section. 

Proof. Since the proof is an adaptation of Chrusciel et al argument, [5J, we 
only sketch the main steps and highlight the point where the non zero force 
is incorporated. We refer the reader to [5J and its references for more details. 

Begin with the expression for the curvature scalar R of the metric g^, 
given by 



l -Re^ = ±Aa + ^(da) 2 + ~A 2 q + ^e~ 2 \pA p , z 



1 

16' 



1 

T 



16 



.4 



z,p; j 



where A is the Laplacian in M 3 and A2 is the 2-dimensional Laplacian 

A 2 g = q, PP + q,zz- 



(16) 



(17) 



Integrate (fTB"j) over M 3 and use the asymptotic flatness condition and the 
definition of the mass m (see |10j) : 



Aad s x 



lim 

r— >oo 



d r ar 2 dS = -8 



dB r 



to obtain 
27rm = 



16 



(da) 2 + -A 2 q + -Re^ + -p\-^\ p A^ - A 



16' 



l z,p) 



d x. 



(19) 

Use inequality (Q to bound the Ricci scalar in the above integral and disre- 
gard the last, non negative term to obtain the following lower bound for the 
mass 



m 



>^f (da) 2 + ^dx* + J- / A 2 qd*x. 
32tt J r3 r] 2 87T J R3 



(20) 



Note that the first integral is exactly the functional Ai used by Dain [10J and 
Chrusciel et al |5j in the proof of the mass-angular momentum inequalities, 



M:=±[ (da) 2 + ^fdx\ 
32tt J r3 i] z 



(21) 



The second integral in the right hand side of (120]) is zero for regular solu- 
tions (i.e, when q = at the axis). But when the solution contains conical 
singularities it gives instead 



„ N-l 

/ A 2 qd 3 x = 2-k QiL 

3 



(22) 
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where Lj is the Newtonian length of Tj and q^ is the constant value of q on 
it. With this result, we arrive at the inequality 

1 N ^ 

i=i 

In [5] (see proposition 2.1) the authors prove that for any set of aligned 
punctures and of axis values Ui := Co>|rv, there exists a solution {a, u) of the 
variational equations associated with the functional Ai, with finite values of 
M. and appropriate asymptotic behavior near each puncture. So, denote by 
f(Ji, --Jn) the numerical value of M of the harmonic map, (a, u) from R 3 \T 
to the two dimensional hyperbolic space. Then we obtain 

1 N ^ 

m>f{J 1 ...J N ) + -^2g i L i . (24) 

i=l 

Finally expressing the constants qi in (124"|) in terms of the forces (see equation 
(ITT]) ), the statement of the theorem is proven. □ 



2.2 Bound in terms of area and angular momentum 

In this section we obtain a bound on the forces acting on a black hole in 
terms of its area and angular momentum. This will generalize inequality ([2]) 
to the non-equilibrium case. 

For that purpose, we will follow Dain and Reiris's argument [SJ, using 
minimal stable surfaces. It is remarkable that since this result is quasi-local, 
there is no asymptotically flatness requirement concerning the initial data. 
We only need the existence of stable minimal surfaces located at the axis, 
with aligned angular momenta (in order to maintain axial symmetry). It is 
clear that in order for these surfaces to have a non vanishing Komar angular 
momentum a non trivial topology is needed. One way to think about this is 
imagining that each surface encloses a puncture. 

A surface E is called minimal if its mean curvature x vanishes, and it 
is called stable if it is a local minimum of the area functional ([8]). This 
can be precisely stated in the following form. Consider a flux of surfaces 
F t : R x S 2 -> S parametrized by t G R such that F| t=0 (S' 2 ) = E and 
F a \ t=0 = an a where dot denotes derivative with respect to t, n a is the unit 
normal to E and a is in principle an arbitrary function on E. Then the 
stability condition can be writen as 

A\ t=Q = / axdSy, (25) 
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where 



x = _A 7 a -(R-R + X abX a >, (26) 



Xab is the second fundamental form of S, R is the Ricci scalar associated to 
the 3-metric g a b (see equation and A 7 and R 1 are the Laplace operator 
and the Ricci scalar associated with the 2-metric jab- 

As in the previous section, we allow the function q in the 3-metric to 
have non zero values at the bounded components Tj. This leads us to the 
following result, which is an extension of theorem 1 in [8j. 

Theorem 2.2. Consider an axisymmetric, vacuum and maximal initial data, 
with a non negative cosmological constant and N > 1 orientable closed stable 
minimal axially symmetric surfaces Then 

v /l + 4J- l _ lV / l + 4J- l >^4^ i = l,..N (27) 

where Ji and are the angular momentum and area of Ej respectively. 

Before giving the proof of theorem 12. 2\ we want to remark that the equal- 
ity in (127|) is achieved by a solution describing two counter-rotating identical 
extreme Kerr black holes separated by a massless strut given by Manko et 
al, (2TJ. The explicit values for the separation distance, ratio A/8tcJ and the 
force in the extreme case as 

L 2 a 2 + m 2 A 2\J\ ^_ (a 2 - m 2 ) 2 

a 2 — m 2 ' 87r|J| a 2 + m 2 ' 16 J 2 

where J = am and J is the angular momentum of each black hole, and m its 
mass parameter (also, 2m is the total ADM mass). We note that a 2 > m 2 
holds for the extreme counter-rotating constituents, unlike in the case of a 
single Kerr black hole, for which a 2 = m 2 . It is also clear from the above 
expressions that as m — > a~, then the separation distance becomes arbitrarily 
large, and the mutual force becomes arbitrarily small. 

In case that the force is positive, which is the expected outcome, the 
quantity A/87r|J| can be less than unity, a situation that can not happen 
for regular black holes. Nevertheless, a type of variational characterization 
of extreme Kerr throat is still valid in the presence of a conical singularity. 
This is due to the fact that in the variational principle used to derive the 
above inequality (see [I]), the function q plays no rol, the relevant functions 
are the conformal factor e a and the twist potential u. 

The inequality ([2]) has been recently extended in [13J, [14| to include elec- 
tric and magnetic charges and to apply also to stable marginally outermost 
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trapped surfaces satisfying the dominant energy condition. To be more pre- 
cise, consider a surface E embedded in a spacetime, with null normal vectors 
£ a and k a such that £ a k a = — 1 and £ a is outward pointing. Let 9^ be the 
expansion associated with the null normal i a . Then, E is a stable marginally 
outermost trapped surface if 9^ = and if there exists an outgoing vector 
X a = j£ a — ipk a with 7 > and ij) > such that 8x0™' > 0. Here 5 X denotes 
the deformation operator on E that controls the infinitesimal variations of 
geometric objects defined on E under an infiniteimal deformation of the sur- 
face along the vector X a . By following the same lines as in theorem (12. 2p 
based on the proof of j 13 j . |14] one obtains the following corollary 

Corollary 2.3. Consider an axisymmetric, stable marginally outermost trapped 
surface Ej ; with a non negative cosmological constant, satisfying the domi- 
nant energy condition. Then 

(29) 

where Ji, Qie, QiB o,nd Ai are the angular momentum, electric charge, mag- 
netic charge and area o/Ej respectively. T{ and J r i _i are the forces acting on 
Ej along the axis componentes Tj and IV_i. 

In the case of zero angular momenta, the inequality is saturated for the 
Majumdar-Papapetrou solution with Ji = QiB — 0, Ai — 4nQ 2 E and van- 
ishing forces. The same thing does not seem plausible for Kerr black holes, 
since numerical results [9j suggest that a configuration with two Kerr black 
holes having a = m leads to a positive force. 

Proof of theorem \2. 6 A We follow the lines of [8] and refer the reader to this ar- 
ticle for more details. Choose a coordinate system such that the determinant 
of the induced metric jab on Ej is 

y/det{i) = e Cl sin 9 (30) 

where c« is a constant. Then, the metric can be written as 

lab = e ai (e 2m (dp 2 + dz 2 ) + p 2 sin 2 9d<f) 2 ) (31) 

with <Tj + qi = Ci and 

(o~i + qi)\e=o = {?i + qi)\e=n = c% (32) 

Now multiply x, given by (1261) by a, use the constraint to express the 
Ricci scalar R in terms of K ao , integrate over Ej and use condition x = 0, 
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the stability condition (I25p and the bound Then we find, disregarding 
the non negative terms 



/ (\Da\ 2 + -R y a 2 )e c dS > [ -^e a ^a 2 e c ^dS, 
J 2 J 4 rj 2 



(33) 



where a prime denotes derivative with respect to 9. Now, using that 
sin 6 



i? 7 = __ (2q' i cos 9 + sin Qo\q\ + 2 sin - (sin Oa'J) (34) 



and choosing 

a = e c ^ /2 (35) 

we obtain 

47r ( c . + i) _ 27 r(g i ( 7 r) + ft (0)) " /(^ + )^ > | (36) 

which can be written as 

4tt( Ci + 1) - 27r( gi (7r) + ft (0)) > (37) 

where A^i is the functional introduced by [1] and used also by [8] and [17] . 

Mi := 7^ / 4o-i + of + (38) 

ZTT J Tj 

Note that in going from f )33|) to f )36|) is where we have used the non regularity 
of the metric at the axis. Now, in [lj it has been proven that for each Mi, 
the following bonds holds 

Mi > 8(ln(2|Ji|) + 1), (39) 

where Jj is the angular momentum of Ej. Using the fact that = 4ire Ci we 
find 

3,-(0)+9,-(7r) 

At > 8n\Ji\e * (40) 

Finally, using the expression for Ti in terms of qi we complete the proof of 
the theorem. 

□ 
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